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ABSTRACT 

Aims. Abundances of 18 chemical elements have been derived for 11 A (normal and Am) and 11 F dwarfs members of the Coma 
Berenices open cluster in order to set constraints on evolutionary models including transport processes (radiative and turbulent dilfu- 
sion) calculated with the Montreal code. 

Methods. A spectral synthesis iterative procedure has been applied to derive the abundances from selected high quality lines in high 
resolution high signal-to-noise echelle spectra obtained with ELODIE at the Observatoire de Haute Provence. 

Results. The chemical pattern found for the A and F dwarfs in Coma Berenices is reminiscent of that found in the Hyades and the 
UMa moving group. In graphs representing the abundances [X/H] versus the effective temperature, the A stars often display abun- 
dances much more scattered around their mean values than the F stars do. Large star-to-star variations are detected for A stars in their 
abundances of C, O, Na, Sc, Ti, Mn, Fe, Ni, Sr, Y, Zr and Ba which we interpret as evidence of transport processes competing with 
radiative diffusion. 

The abundances of Mn, Ni, Sr and Ba are strongly correlated with that of iron for A and Am stars. In contrast the ratios [C/Fe] and 
[0/Fe] appear to be anticorrelated with [Fe/H] as found earlier for field A dwarfs. All Am stars in Coma Berenices are deficient in 
C and O and overabundant in elements heavier than Fe but not all are deficient in calcium and/or scandium. The F stars have solar 
abundances for almost all elements except for Mg, Si, V and Ba. 

The derived abundances patterns, [X/H] versus atomic number, for the slow rotator HD 108642 (A2m) and the moderately fast rotator 
HD106887 (A4m) were compared to the predictions of self consistent evolutionary model codes including radiative and different 
amounts of turbulent diffusion. None of the models reproduces entirely the overall shape of the abundance pattern. 
Conclusions. While part of the discrepancies between derived and predicted abundances may be accounted for by non-LTE effects, 
the inclusion of competing processes such as rotational mixing in the radiative zones of these stars seems necessary to improve the 
agreement between observed and predicted abundance patterns. 

Key words, stars: abundances - stars: main sequence - stars: rotation - Diffusion - Galaxy: open clusters and associations: individual: 
Coma Berenices 

1 . Introduction determinations of 18 chemical elements in 11 A and IIP dwarfs 

in the Coma Berenices open cluster (age about 447 Myrs). The 
Abundance determinations of A and F dwarfs in open clusters spectroscopy presented here is part of an ongoing observational 
and moving groups of known properties aim at elucidating program of A and F dwarfs in galactic open clusters of various 
the mecanisms of mixing at play in the interiors of these ^ges whose aims are twofold. First, we wish to improve our 
main-sequence stars. In two previous papers, abundance knowledge of the chemical composition of A dwarfs (normal A 
determinations were presented for 11 chemical elements in stars and Am stars) which is still poor in particular for normal 
the Hyades (age about 787 Myrs, [Varenne & Monier 1999|l a stars. Second, we intend to use the derived abundances to 
andjnjhe Ursa Major moving group (age about 500 Myr, set constraints on self-consistent evolutionary models of these 
|Monier2005). Specifically, 19 A and 29 F dwarfs members of objects including various particle transport processes. Indeed 
the Hyades and 12 A and 10 F dwarf bona-fide and probable stars in open clusters originate from the same interstellar ma- 
members of the Ursa Major moving group were analysed. Serial (ie they have the same age and the same initial chemical 
The motivation of this study is to report on new abundance composition) and as such are very useful to test the predictions 
of evolutionary models. 

Send offprint requests to: M. Gebran F^w Studies have addressed the chemical composition of the 

* Based on observations at the Observatoire de Haute-Provence A and F dwarfs in the Coma Berenices open cluster High 
(France). quality high resolution spectroscopy is feasible for the brightest 

** Present aflfliation: Laboratoire Universitaire d' Astrophysique de members because of the proximity (d ^ 80 -90pc) of the cluster. 
Nice, UMR 6525, Universite de Nice - Sophia Antipolis, Pare Valrose, Abundance analyses of the Coma Berenices dwarfs have mostly 
06108 Nice Cedex 2, France 
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focused on F and G dwarfs whose low apparent rotational 
velocities facilitate the abundance determination. For concision, 
we have collected previous abundance determinations of the A, 
F and G dwarfs in Coma Berenices in Table [1] where the biblio- 
graphical references, the number of stars analysed, their spectral 
types and the investigated chemical elements are collected. 
Lithium abundances have been determined for several F (and G) 
dwarfs in Coma Berenices by Boesgaard (1987) Jeffries (1999)| 



and [Soderblom et al. (1990)| Carb on and iron abundances 



were derived for 14 F dwarfs by Friel & Boesgaard (1992) 
Carbon, oxygen, silicon, calcium, iron, barium, magnesium 
scandium, chromium, nickel, lithium, aluminium, sulfur and 
europium abundances have been determined for only a few 
normal A and Am stars (fewer than 7 stars) by Savanov (1996), 
Hui-Bon-Hoa et al. (1997), H ui-Bon-Hoa & Alecian (1998) 
and Burkhart & Coupry (2000)| The abundance determinations 
for A stars focuse mainly on the chemically peculiar Am stars. 
[Savanov (1996)| found significant star-to-star differences in 
abundances among the Am stars of Coma Ber for a given 
chemical element. In contrast, little attention has been paid to 
the "normal" A stars of the cluster The chemical composition 
of normal A dwarfs (field and cluster stars) remains generally 
poorly known as too few objects have been analysed so far, 
mainly because of their high rotational velocities. Significant 
abundance differences have been found among the few field 
"normal" A stars analysed so far (Holweger et al. 1986 



ILambertet al."7986l Lemke 1998,1990; Hill & Landstreet 1993; 
Hill 1995; Rentzsch-Holm 1997; Varenne 1999). Varenne & 
Monier (1999) also found significant star-to-star variations of 
the abundances of O, Na, Ni, Y and Ba for the normal A and the 
Am stars in the Hyades whereas the F dwarfs display much less 
scatter. Similarly, Monier (2005) found star-to-star variations in 
[Fe/H], [Ni/H] and [Si/H] much larger for the A dwarfs than for 
the F dwarfs in the Ursa Major group. 

The main thrust of this paper is to report on the abundances 
of 18 chemical elements (C, O, Na, Mg, Si, Ca, Sc, Ti, V, 
Cr, Mn, Fe, Co, Ni, Sr, Y, Zr and Ba) for F and A dwarfs in 
Coma Berenices. These abundances have been compared to 
published predictions of recent self consistent models including 
transport processes at ages close to that of Coma Berenices 
(Turcotte et al. 1998, Richer et al. 2000, Richard et al. 2002) 
or to new models calculated with the Montreal code (Richard 
et al. 2002). The abundances were derived by synthesizing 
carefully selected lines in high quality high resolution spectra 
of 1 1 A and Am stars and 1 1 F stars of the cluster The selection 
of the sample of stars observed, the observations and the data 
reduction are described in Section 2. The determination of the 
fundamental parameters, the construction of the linelist and the 
spectrum synthesis are discussed in Section 3. The behaviour 
of the abundances of individual chemical elements in the A 
and F dwarfs of Coma Berenices versus effective temperature 
and the abundance of iron are described in Section 4. The 
astrophysical implications of our findings and the detailed 
comparison of the found abundance patterns for 2 Am stars with 
recent self-consistent evolutionary models including transport 
processes are presented in Section 5. Conclusions are given in 
Section 6. 



2. Program stars, observations and data reduction 

Our observing sample consists of all A and F stars members of 
the Coma Berenices cluster brighter than V - 8.6 mag retrieved 
from Trumpler's (1938) list. This magnitude corresponds to the 



latest F dwarfs (F8V) which required the longest achievable ex- 
posure times (about 1 hour 15 minutes). Eleven A ("normal" and 
Am stars) stars and 1 1 F stars members of the Coma Berenices 
cluster were observed using two spectrographs at Observatoire 
de Haute Provence (OHP): the ELODIE echelle spectrograph 
and the monoorder AURELIE spectrograph. AURELIE is a 
monoorder spectrograph (Gillet et al. 1994) placed at the coude 
focus of the 152 cm telescope. AURELIE is free of scattered 
light and has been used to calibrate putative effects of scat- 
tered light on the line profiles in the blue region obtained with 
ELODIE. ELODIE is a fiber-fed cross-dispersed echelle spec- 
trograph attached to 193 cm at OHP (Baranne et al. 1996). It 
records in a single exposure a spectrum extending from 3850 A 
to 6811 A at a resolving power of about 42000 on a relatively 
small CCD (1024X1024). The observing dates, exposure times 
and Signal to Noise ratios achieved for each star are collected in 
Table [3] the first part being dedicated to ELODIE spectra while 
the second describes the AURELIE spectra obtained in three 70 
A wide regions centered around 4505 A, 5080 A, 5530 A and 
6160 A. 

The fundamental data of these stars are collected in Table 
|2] The Trumpler and Henry Draper identifications appear in 
columns 1 and 2, the spectral types retrieved from SIMBAD in 
column 3, the apparent magnitudes in column 4. Columns 5 to 8 
display the effective temperatures and surface gravities adopted 
for the analysis and the rotational velocities and microturbulent 
velocities derived in our analysis (see section [J.2.1t . Comments 
about binarity and pulsation appear in the last column. Note that 
2 stars, HD 107655 and HD 106999, may not be members of 
the cluster (Bounatiro 1993). Although we derived abundances 
for these stars too, their data do not appear in the figures. There 
are no very rapid rotators in this cluster, the apparent rotational 
velocities range from 9 km-s"' to 102 km-s"'. 
Inspection of the Catalogue of Double and Multiple stars 
(CCDM, Dommanget & Nys 1995) reveals that only one star is a 
binary system. HD 106887 is the primary star in a double system 
and has a much fainter companion (V - 9.8) located nearby at 
about 8.6 arcsec, whose spectral type is unknown. It is unlikely 
that the light of the companion may have contributed signifi- 
cantly to the observed spectra of HD 106887. 
With AURELIE, four spectral regions centered on /14505 A (re- 
gion 1), /16160A (region 2), /15080A (region 3) and ^5530A (re- 
gion 4) have been observed as they include several lines having 
accurate oscillator strengths for the chemical elements we study. 
Region 1 contains the Mg 2 triplet at /1448 1 A and several clean 
unblended and well separated Fe 2 lines. Regions 2 to 4 are those 
observed by Edvardsson et al. (1993) in their spectroscopic sur- 
vey of F dwarfs in the galactic disk. For the brightest stars (V < 
6), grating 5 working at its second order (for region 1) and grat- 
ing 1 (for region 2) yielded resolving powers equal to 60000 
and 65000 respectively. For the fainter stars (V > 6), grating 
7 was used, resulting in resolving powers of 36000 and 28000 
(for regions 1 and 2 respectively). In region 3, all stars have 
been observed with grating 7 (R = 31000) except for Procyon. 
According to the V magnitude of the stars and weather condi- 
tions, exposure times between lh30 and 3 hours were necessary 
to secure signal to noise ratios around 200. The IRAF software 
was used to reduce all spectra following the standard procedure 
(mean offset removal, division by the mean flat field, wavelength 
calibration and continuum normalization). 
At the end of an observing night, the ELODIE spectra 
are provided to the observer fully reduced by the INTER- 
TACOS (INTERpreter for the Treatment, the Analysis and the 
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Table 1. Previous abundance determinations for the Coma Ber A, F and G dwarfs. 



Reference 


Stars studied 


Chemical Elements 


Savanov (1996) 


13 A- Am and F-Fm 


C,0,Si,Ca,Fe,Ba 


Hui-Bon-Hoa et al. (1997) 


2 A-Am 


Mg,Ca,Sc,Cr,Fe,Ni 


Hui-Bon-Hoa & Alecian (1998) 


4 A-Am 


Mg,Ca,Sc,Cr,Fe,Ni 


Burkhart & Coupry (2000) 


7 A-Am 


Li,Al,Si,S,Fe,Ni,Eu 


Boesgaard (1987) 


22 A and F 


Li 


Friel & Boesgaard (1992) 


14 F 


Fe,C 


Jeffries (1999) 


15 F,G,K 


Li 


Cayreletal. (1988) 


4G 


Fe 


Soderblom et al. (1990) 


28 G 


Li 



Table 2. Data on the programme stars. Spectral type are taken from SIMBAD and WEBDA online database. Tgff and \ogg are 
those determined by UVBYBETA code. Vg sin / and ^, are determined as explained in sect. 13.2.11 References (a) and (b) are 
[Bounatiro (1993)| and Rodriguez et al. (1994) respectively. 



TR 


HD 


Type 


m,, 


(K) 




sin I 
km-s"' 


km-s ' 


Remarques 


19 


HD106103 


F5V 


8.09 


6707 


4.45 


19.7 


1.4 






HD 106293 


F5V 


8.09 


6545 


4.34 


47 


1.4 




36 


HD 106691 


F5IV 


8.08 


6713 


4.43 


37 


1.6 




49 


HD 106946 


F2V 


7.87 


6892 


4.30 


62 


1.9 




86 


HD107611 


F6V 


8.50 


6491 


4.57 


22 


1.4 






HD109530 


F2V 


7.30 


6497 


3.85 


67 


2.1 




101 


HD107877 


F6 


8.35 


6598 


4.54 


29.5 


1.25 




114 


HD108154 


F5 


8.56 


6497 


4.54 


19 


1.15 




118 


HD108226 


F5 


8.34 


6530 


4.45 


19 


1.05 




162 


HD 108976 


F6 V 


8.54 


6413 


4.49 


20 


1.1 






HD 109069 


FO V 


7.55 


6864 


4.06 


89 


2.2 





107 


HD 107966 


A3V/A31V 


5.18 


8541 


3.82 


51 


2.9 




130 


HD108382 


A4V/A31V 


4.96 


8317 


3.92 


75.5 


3.2 




47 


HD106887 


A4m 


5.71 


8291 


4.20 


82 


3.8 




86 


HD107655 


AOV 


6.18 


9675 


4.10 


45 


2.0 


^ Coma (a) 


62 


HD107168 


A8m/kA5hA5mFO 


6.24 


8283 


4.20 


14.3 


4.0 




183 


HD109307 


A4Vm/A3IV-V 


6.26 


8396 


4.10 


14.5 


3.3 




144 


HD 108642 


A2m/kA2hA7mA7 


6.54 


8079 


4.06 


9.2 


3.7 




68 


HD 107276 


Am/Ka5mA7 


6.63 


8000 


4.00 


102 


2.8 




139 


HD 108486 


AmkA3hA5mA7 


6.67 


8148 


4.11 


37 


3.0 




52 


HD 106999 


Am 


7.46 


8148 


4.09 


44 


3.0 


i Coma (a) 


82 


HD107513 


Am/kA7hF0mF0 


7.38 


7279 


4.02 


62 


3.0 


d Scuti (b) 



correlation of Spectra) pipeline developed by D. Queloz and L. 
Weber (Baranne et al. 1996). The reduction is actually included 
in the spectrograph data flow and provides fully flat-fielded and 
wavelength calibrated spectra directly at the end of each expo- 
sure. However, we have chosen to perform our own reduction of 
the ELODIE spectra. Indeed, ELODIE was primarily designed 
to provide accurate radial velocity measurements. One source 
of concern when deriving abundances is to properly correct the 
spectrum for scattered light, especially in the blue region. The 
background in a stellar ELODIE exposure can be estimated by 
measuring the flux in the inter-orders. In INTER-TACOS, this 
background is removed using a two dimensional polynomial fit 
with a typical 5 % error which peaks in the middle of the orders 
(see fig. 11 in Baranne et al. 1996). Erspamer & North (2002) 
have devised a reduction procedure using a set of IRAF func- 
tions which they apply to the raw image in order to provide 
an improved correction for scattered light compared to INTER- 
TACOS. We have also chosen to perform our own reduction of 



the ELODIE spectra using IRAF (Image Reduction and Analysis 
Facility, Tody 1993) routines. Although it follows Erspamer & 
North's (2002) procedure, our reduction slightly differs from 
their method. The steps are as follows: 

1- Averaging the several offsets and flat-fields taken throughout 
the night, using zerocombine and flatcombine. 

2- Removal of the mean offset from all images and bad pixels 
correction using ccdproc . 

3- Finding and centering the orders using the flat-field image 
using apfind and apcenter. 

4- Removal of the scattered light (apscatter). Scattered light 
fills in the line profiles, making the lines shallower and 
thus leading to underabundances if not taken into account. 
The scattered light is estimated (as explained in sect 3.2 of 
Erspamer & North 2002) and substracted from the original 
image. An example of the effect of this improved removal of 
the scattered light is shown in Figure |2] where the corrected 
profiles of the Fell fines at 4520.224 A and 4522.634 A are 



4 



M. Gebran, R. Monier and O. Richard: Chemical composition of A and F dwarf members of the Coma Berenices open cluster 



compared to the INTER-TACOS profiles. In this case, ignor- 
ing the scattered light would lead to underestimation of the 
iron abundance deduced from these 2 lines by 0.07 and 0.08 
dex respectively. The effect is more pronounced, about 0.14 
dex, for the Mgll triplet at 4481 A. 

5- Extraction of the images using apsum. The averaged flat- 
field is used to determine the shape of the orders which is 
used later as a reference for the extraction of the images. 
The extraction method uses Home's (1986) algorithm. 

6- Calibration of the thorium image using apsum, ecidentify 
and ecreidentify. 

7- Division by the extracted flat-field using sarith. 

8- Wavelength calibration of the spectra using the thorium spec- 
tra using dispcor. 

9- Normalization to the continuum using continuum. To ensure 
we correctly located regions free of lines (when available) in 
each order, we have computed synthetic spectra using the 
code SYNSPEC48 (Hubeny & Lanz 1992) assuming a solar 
metallicity for the various temperatures and surface gravities 
of our stars. The spectrum was then rectified to this local 
continuum. 

10- Merging the 67 normalized orders using scombine. The last 
three orders do not overlap and thus could not be merged. 
For the abundance analysis, we have discarded lines located 
in the overlapping region of two successive orders. 

11- Radial velocity determination using fxcor: the final merged 
spectrum is cross-correlated with differents masks of spectral 
types AOV, A5V, A9V and F5V to derive the radial velocity. 
The merged spectrum is then corrected for the radial velocity 
found. 

Figure [3] compares the spectral order 21, centered around A45QQ 
A, for 4 A and 4 F stars reduced in this manner. The effect of 
increasing stellar rotation on line profiles is conspicuous. The 
last plot of Figure [3] displays a typical agreement between the 
observed spectrum of HD 107655 (AOV) and the synthetic spec- 
trum (computed as explained in sect. O that provides the best 
fit. 



3. Abundance analysis 

The abundances of 18 chemical elements have been derived by 
iteratively adjusting synthetic spectra to the observed normalized 
spectra and minimizing the chisquare of the models to the ob- 
servations. Spectrum synthesis is the most appropriate method 
as our stars have apparent rotational velocities ranging from a 
9 to 102 km-s"'. Specifically, synthetic spectra were computed 
assuming LTE using Takeda's (1995) iterative procedure and 
double-checked using Hubeny & Lanz's (1992) SYNSPEC48 
code. This version of SYNSPEC calculates lines for elements 
heavier than Zn up to Z=99. 

3.1. Atmospheric parameters arid model atmospheres 



The effective 
ties (logg) of 



Napiwotzki et al. (1993) 



temperatures (Teg) and surface gravi- 
the stars have been determined using the 
UVBYBETA caHbration of the 



Stromgren photomery indices uvby in terms of Teff and logg. 
The found effective temperatures and surface gravities are 
collected in Table |2] The errors on Teff and logg are estimated to 
be + 125 K and ± 0.2 dex respectively. Model atmospheres were 
then calculated using Kurucz's ATLAS9 code (Kurucz, 1979), 
assuming a plane parallel geometry, a gas in hydrostatic and ra- 
diative equilibrium and LTE. The ATLAS9 model atmospheres 



BOOO 



Fig. 1. Hertzsprung-Russel diagram of the observed stars in 
Coma Berenices (visual magnitude versus T^ff). The two stars 
flagged with question marks are probably not members of the 
cluster. 




Fell (0.07 dex) 



Fell (0.08 dex) 



.... T.4C0S reducLio 
IRAr rcdiiotion 



4516 4518 45S0 



15^4 4526 4528 4530 



Fig. 2. Effect of a proper removal of scattered light. The cor- 
rected profiles (IRAF reduction, thick lines) are deeper than the 
INTER-TACOS ones (dashed lines). The abundances deduced 
from the corrected profiles are about 0.08 dex larger. 



contain 64 layers with a regular increase in log tross = 0. 125 and 
they have been calculated assuming Grevesse & Sauval (1998) 
solar chemical composition. This ATLAS9 version uses the new 
opacity distribution function (ODFs) (ICastelli & Kurucz 2003l l 
calculated for this solar chemical composition. The line opacity 
calculation uses the 58 million linelist compiled by Kurucz 
(Kurucz 1992a, 1992b). Convection in ATLAS9 relies on the 
mixing length theory (MLT). Specifically, we have adopted 
Smalley's prescriptions ( jSmalley 2004| l for the values of the 
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o I 1 I 1 I 1 I 1 I 1 I 1 I 

4480 4490 4500 45 lO 4520 4530 



1.2 



0.8 



0.6 



Observed spectrum 
Synthetic spectrum 




HD107655 (AOV) 

V sini— -45 km/s _ 



_L 



4480 



4490 



4500 



45 lO 



4520 



Fig. 3. Selected observed spectra of A (top) and F (middle) dwarf stars members of Coma Berenices open cluster. Spectra are 
arbitrarily shifted vertically by 0.5, 1 and 1.5 unit of normalized flux. The smearing out of the spectra by rotation is noticeable. 
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Table 3. Observing log of the programme stars of Coma Berenices. The first table describes the ELODIE observations and the 
second one the AURELIE observations 



HD 


spectral 


Mv 


exposure 


S/N 


Date 




type 


mag 


time (s) 






106103 


F5V 


8.09 


4500 


144 


04/10/04 


106293 


F5V 


8.09 


4500 


132 


04/10/04 


106691 


F5IV 


8.08 


4500 


147 


04/10/04 


106946 


F2V 


7.87 


4500 


180 


04/11/04 


107611 


F6V 


8.50 


4500 


129 


04/11/04 


109530 


F2V 


7.30 


4500 


185 


04/12/04 


106887 


A4m 


5.71 


3600 


283 


04/07/04 


106999 


Am 


7.46 


4500 


177 


04/10/04 


107131 


A6IV-V 


6.44 


4500 


213 


04/09/04 


107168 


A8m 


6.24 


4500 


215 


04/07/04 


107276 


Am 


6.63 


4500 


152 


04/09/04 


107513 


Am 


7.38 


4500 


184 


04/10/04 


107655 


AOV 


6.18 


4500 


254 


04/07/04 


107966 


A3V 


5.18 


3600 


336 


04/07/04 


108382 


A4V 


4.96 


3600 


337 


04/07/04 


108486 


Am 


6.67 


4500 


121 


04/09/04 


108642 


A2m 


6.54 


4500 


170 


04/09/04 


108651 


AOp 


6.65 


4500 


450 


04/11/04 


109307 


A4Vm 


6.26 


3600 


114 


04/08/04 



HD 


spectral 


Mv 


grating 5: /lc=4505 A 


Date 


grating 7: /lj.=6160 A 


.{,=5080 A 


.{,=5530 A 




type 


mag 


exp.time(s) 




exp.time(mn) 


exp.time(mn) 


exp.time(mn) 


107168 


A8m 


6.24 


800 


03/14/04 


75 


104 


120 (grating 1) 


107513 


Am 


7.38 


2000 


03/15/04 


120 


120 


120 


106103 


F5V 


8.09 


4500 


03/15/04 


60 






107655 


AOV 


6.18 


800 


03/15/04 


135 


85 


90 


106887 


A4m 


5.71 


600 


03/15/04 






60 


107131 


A6IV-V 


6.44 


800 


03/15/04 


110 






107276 


Am 


6.63 


1000 


03/15/04 


90 


120 


75 


107966 


A3V 


5.18 


500 


03/15/04 


45 


105 


40 


108382 


A4V 


4.96 


500 


03/15/04 


60 


74 




108486 


Am 


6.67 


1000 


03/15/04 




180 


90 


108642 


A2m 


6.54 


1000 


03/15/04 


90 


120 




108651 


AOp 


6.65 


1000 


03/15/04 


105 


150 




109307 


A4Vm 


6.26 


800 


03/15/04 


120 


120 (grating 1) 


80 


106293 


F5V 


8.09 


4500 


03/16/04 


130 






106691 


F5IV 


8.08 


4500 


03/16/04 


120 


150 


120 


109069 


FOV 


7.55 


3000 


03/16/04 


100 


150 


110 


107877 


F6 


8.35 


4500 


03/16/04 


180 


165 




107611 


F6V 


8.50 


4500 


03/16/04 


90 






108154 


F5 


8.56 


4500 


03/17/04 


180 




165 


108226 


F5 


8.34 


4500 


03/17/04 


90 


180 




108976 


F6V 


8.54 


4500 


03/17/04 


150 




140 


109530 


F2V 


7.30 


3000 


03/17/04 


100 


120 


110 



ratio of the mixing length to the pressure scale height (a - ^) 
and also for the microturbulent velocity (constant with depth). 

3.1.1. The linelist 

The linelist was constructed from Kurucz's gfall.dad list, from 
which we selected lines between 3000 and 7000 A. However the 
data in this list have been modified and complemented in differ- 
ent manners. For the lines which we expected to contribute sig- 
nificantly to the absorption, we have carefully checked the wave- 
lengths, lower excitation potential, oscillator strength and damp- 
ing constants (radiative. Stark and Van der Waals) in gfall.dat 
against more accurate and/or more recent critically evalu- 
ated laboratory determinations when available. Specifically, two 
atomic databases were searched for improved values of these 

1 http://kurucz.harvard.edu/LlNELlSTS/GFALL/ 



parameters and their uncertainties: the VALEQ database and the 
NIS10 database. We then modified their values in the original 
linelist accordingly. Damping constants not available in linelists 
are calculated in SYNSPEC48 using approximations. We have 
also excluded lines in the overlaping regions of two succes- 
sive orders, as we felt that the observed line profiles may not 
be reliable there. The final linelist contains 270 transitions for 
18 elements which we believe are reliable enough to derive the 
abundances. Most of the lines studied here are weak lines which 
are formed deep in the atmosphere where LTE should prevail. 
They are well suited for abundance determinations. The final 
linelist appears in Table |8] where, for each element, the wave- 
length, adopted oscillator strength, its accuracy (when available) 
and original bibliographical reference are given. We have also 
included data for hyperfine splitting for the selected transitions 

" http ://ams . astro .univie. ac . at/vald/, |Kupka et al. 1999] 
' http://physics.nist.gov/cgi-bin/AtData/lines-form 
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when relevant, in particular for Mn II (these were retrieved from 
the linelist gfhyperall.dafl). However the moderate spectral res- 
olution of the spectra and smearing out of spectra by stellar ro- 
tation clearly prevent us from detecting signatures of hyperfine 
splitting and isotopic shifts in our spectra. 

3.2. Synthetic spectra 

Two codes have been used to derive abundances: Takeda's code 
dTakeda 1995 | and SYNSPEC (Hube ny & Lanz 1992 1) to check 
the abundances produced by Takeda's iterative procedure. We 
briefly review the assumptions of these codes in section 13.2.11 
and section [J.2.2| 

3.2.1. Takeda's procedure 

Takeda's procedure iteratively minimizes the dispersion cr be- 
tween the normalized synthetic spectrum and the observed one 
defined as: 

^2 _ {y,-n,-Cf 

where is the number of wavelength points, y, the logarith- 
mic of the observed spectra (f^.) and 77, the logarithmic of the 
synthetic spectra iF,\.). C is an offset constant reflecting a pos- 
sible difference of units between and (it should be very 
close to zero when working with normalized fluxes). The syn- 
thetic flux 77, in a line is a priori a function of several physical 

variables (xi , xk) which represent the unknowns one might be 

looking for: abundances of individual elements, projected rota- 
tional velocity and microturbulent velocity, oscillator strengths, 
damping constants. The dispersion cr^ is thus a function of K-hl 
variables (including C). Minimizing cr^ thus requires that its 
partial derivatives with respect to these K-i-1 variables be zero, 
which leads to a set of K-i-1 non-linear equations with K+1 un- 
knowns. Takeda (1995) solves this problem numerically using 
a Newton-Raphson iterative technique (linearization method). 
Iterations are repeated until convergence ie. when variations in 
the x^s become sufficiently small, typically less than 10""'. 
Takeda's code consists of two routines: the first is a modified 
version of Kurucz's Width9 code dKurucz 1992a l for computing 
the opacity data which needs a Kurucz's model atmosphere as 
input. Once the opacities are computed, the second routine com- 
putes the emergent flux and minimizes the dispersion between 
the synthetic and the observed spectra. 

In this study, we kept the number of free variables at 3. At each 
iteration, the code outputs 3 parameters: the abundance(s) of the 
studied element(s) (log e), the microturbulent velocity (^,) con- 
stant whith depth and the projected rotational velocity (vv sin /), 
the oscillator strengths and damping constants being kept fixed. 
We usually synthesized unblended lines of one chemical element 
only, although the procedure allows us to simultaneously synthe- 
size lines of different elements. 

Before tackling the abundance of each individual chemical el- 
ement, we derived the rotational and microturbulent velocities 
for each star as follows. Allowing small variations around solar 
abundances of Mg and Fe, we iteratively fitted the unblended 
line profile of the Mgll triplet at 4480 A and a set of neigh- 
boring unblended weak and moderately strong Fell lines such as 
4491.405 A and 4508.288 A leaving ^, and vv sin / as free pai-am- 
eters. Even in the fastest rotators, the continuum is well seen in 
this spectral region where there are only a few regularly spaced 

http://kurucz.harvard.edu/LINELISTS/GFHYPERALL/ 
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Fig. 4. Iterative adjustment of synthetic profiles of the Fell line at 
4508.288 A for the A3V star HD107966 with Takeda's program. 
Convergence was achieved properly at the fifth iteration. 

lines. The weakest Fell lines are sensitive to rotational veloc- 
ity and not to microturbulent velocity, the moderately strong 
Fell lines respond to microturbulent velocity changes. The Mgll 
triplet, commonly used to derive stellar rotational velocities, ap- 
pears to be sensitive to both ^, and vv sin /. Each line of Fe and 
Mg yielded a set of values for log e, ^, and sin /, which were 
in good agreement. For the F stars, the found ^, were checked 
against Nissen's (1981) prediction which invoked a polynomial 
fit for as a function of T^ff and log g for the F stars. The agree- 
ment is generally good, the maximum difference we found being 
A^i - 0.32 km-s"'. Morever, the derived VgSini in this analysis 
correlate well with those derived from AURELIE spectra of the 
region around 4500 A published in [Monier & Richard (2004)| 
For the slowest rotating Am stars, curves of growths of the Fe 
II lines were also performed. They led to values of ^, consistent 
with those derived from the fit of the Mg II triplet and Fe II lines 
around 4500 A. Once ^, and sin / were fixed, we then pro- 
ceeded to determine the abundance for each selected unblended 
line of each chemical element. In practice, complete conver- 
gence was reached for each star after up to 10 iterations (five 
in the most favorable cases). 

A final test on Procyon was also performed, assuming a so- 
lar composition ( Steffen 19851) . The iterative fitting of the Mg 
II and Fe II lines around 4500 A yielded an apparent rota- 
tional velocity of 6 km-s"' and a microturbulent velocity of 2.2 
km s ' in good agreement with Steffen's jStefFen 19851 ) values 
(ve sin /=4.5 km-s"' and i;t-2. 1 km-s"' respectively). 
Figure m shows an example of the iterative fitting of the Fell line 
at 4508,288 A for the A3V star HD 107966 (6 iterations). For 
this star, the convergence was already achieved at the fifth itera- 
tion. 

3.2.2. Validation of tine abundances by SYNSPEC 

We have also u sed Version 48 of SYNSPEC 
( [Hubeny & Lanz 1992[ ) to check the abundances produced 
by Takeda's iterative procedure. SYNSPEC48 allows us to 
calculate line profiles of elements up to Z=99. In its LTE mode, 
SYNSPEC needs a model atmosphere and a linelist plus an 
auxiliary file containing non standard flags. For a given star, the 
same ATLAS9 models, the same linelist and the Ve sin / and ^, 
derived using Takeda's (1995) procedure were used to calculate 
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Fig. 5. Influence of the underlying model atmosphere on the de- 
rived abundances for HD107168 (A8m): black dots are abun- 
dances derived using a solar ATLAS9 model, hatched black 
squares using a solar ATLAS 12 model and black triangles are 
abundances derived with a ATLAS 12 model computed for the 
specific chemical composition of this star. 



synthetic spectra. Only the abundance of the studied unblended 
line were left as free parameters. The derived abundances with 
SYNSPEC were found to always agree with those derived from 
Takeda's procedure within the error bars. The abundances listed 
in Tables|4]and|5]are those derived from Takeda's procedure. 
We also checked the influence of the underlying atmo- 
spheric structure on the derived abundances using ATLAS 12 
jKurucz 2005 ). Abundances can be adjusted individually in 
ATLAS12 which employs the Opacity Sampling technique 
for Une opacity. The eff'ect should be noticeable for the Am 
stars whose abundances depart most from the solar values. For 
Am stars we found that the inclusion of an ALTAS12 model 
calculated for the specific chemical composition found with 
ATLAS9 and Takeda's procedure yields new mean abundances 
which diff'er by up to 0.08 dex from those derived with ATLAS9. 
This is less than the eiTor bar on mean abundances. Figure |5] 
compares the abundances of several elements for the A8m 
star HD 107 168 derived using ATLAS9 and ATLAS 12 model 
atmospheres. For the other stars, whose chemical compositions 
depart less from solar, the eff'ect will be even smaller. Since we 
intend to model all stars in a uniform manner and because of 
the much larger computing time needed to run an ATLAS 12 
model atmosphere, we decided to use only ATLAS9 model 
atmospheres. 



4. Results and uncertainties 

Abundances for the 22 stars are collected in Table |4]for A stars 
and in Table |5]for F stars. These abundances are relative to the 
sun ([^] = log(^)* - log(^)0fl They are weighted means of 
the abundances derived from each transition. The calculation 
of their uncertainties and the weighted mean abundances is 
explained in Appendix A. 



' Solar abundances are those from |Grevesse & Sauval (1998)1 



4.1. The general abundance pattern of A, Am and F stars in 
Coma Ber 

Graphs where abundances are displayed against atomic number 
(abundance pattern plots) are particularly appropriate to com- 
pare the behaviour of the A, Am and F stars for different chem- 
ical elements. The abundance patterns for the A stars. Am stars 
and F stars are displayed respectively in Figures |6t,|6j3 and|6};. 
The seven Am stars of the clusters display a characteristic jig- 
saw pattern with much larger excursions from the solar compo- 
sitions than the normal A stars do. This trend had already been 
found by Hui-Bon-Hoa & Alecian (1998) who derived the abun- 
dances of Mg, Ca, Sc, Cr, Fe and Ni in four Am stars of our 
sample (HD107168, HD107966, HD108486 and HD109307). 
Our abundances for these elements agree well with theirs for 
HD 107168, for which all fundamental parameters are simi- 
lar in both studies. However, slight to moderately large dif- 
ferences occur for HD 108486, HD 109307 and HD 107966, for 
which we adopted larger microturbulent velocities (in particu- 
lar for HD 107966), all other parameters being consistent in both 
studies. The 2 normal A stars of our sample, HD107966 and 
HD108382, have very similar abundances, nearly solar, in al- 
most all elements and thus exhibit almost the same abundance 
pattern (Figure [6^). Inspection of Figure |6]3 reveals that gener- 
ally, for almost all chemical elements. Am stars display star-to- 
star variations that can be larger than the typical uncertainty. All 
Am stars are deficient in C and O, but not all are deficient in Ca 
and/or Sc and most have pronounced overabundances of iron- 
peak and heavy elements. The two normal A stars have almost 
solar abundances in most elements. The F stars definitely exhibit 
less scatter than the Am stars. They tend to be mildly overabun- 
dant, in particular in Mg, Si, V and Ba. They have nearly solar 
abundances for C, O, Na, Ti, Fe, Ni, Sr and Y. 

4.2. Behaviour of individual elements 

In this section, we use the found abundances for each element 
to address two issues. First, do the abundances depend on stel- 
lar parameters such as the effective temperature and Vg sin / ? 
Any such correlation could be very valuable for theorists inves- 
tigating the various hydrodynamical mechanisms affecting pho- 
tospheric abundances. Second, how do the abundances of each 
element vary with respect to each other? Given that iron is as- 
trophysically one of the most important elements (since it pro- 
vides a rough estimate of the "metallicity") and given that the 
abundances derived for this element are probably the most re- 
liable, we have examined whether the abundances of individual 
elements correlate with those of iron. Hill (1995) and Lemke 
(1989, 1990) also looked for similar correlations and we will re- 
fer to their findings later. 

We can roughly separate the elements we studied into two 
groups. For most elements (Fe II , Ti II, OI, Cr II, Mg II, Mn 
I, C I, Si II, Ca II and Ni I), we synthesized several lines of 
quality A to D and we can be confident in the abundances we 
derived, in particular for iron, titanium and chromium. For ele- 
ments with many lines, errors in individual oscillator strengths 
should tend to cancel out. For other elements, we have very few 
lines (Sr II), so that errors on oscillator strengths may induce 
scatter in the abundances. For V II, Co I, Y II and Zr II, sev- 
eral lines are available but their accuracy is not necessarily well 
known. The abundances of these elements should be viewed with 
more caution However, our main goal here is to investigate how 
the abundances of individual elements vary with that of iron and 
also to study star-to-star variations of fundamental parameters. 
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Table 4. Abundances relative to hydrogen and to the solar value, [^] - log(^)* - log(^)o for the A stars. The solar values are those 
of |Grevesse & Sauval (1998)] The HD numbers in italics are those for which the uncertainties have been calculated as explained in 
Appendix A. For the others, the quantities labeled as cr are standard deviations. 



HD SpT CI o-c OI o-Q Nal o-a,„ Mgll g-Mg Sill oSj 



HD107966 


A3V/A3IV 


-0.07 


0.04 


0.00 


0.05 


0.52 


0.13 


0.13 


0.09 


0.14 


0.09 


HD108382 


A4V/A3IV 


-0.71 


0.15 


-0.02 


0.06 


0.24 


0.18 


0.13 


0.12 


0.18 


0.19 


HD106887 


A4m 


-0.56 


0.15 


-0.43 


0.06 


0.54 


0.08 


0.17 


0.13 


0.23 


0.20 


HD107655 


AOV 


-0.75 


0.19 


-0.58 


0.06 


0.55 


0.21 


-0.12 


0.18 


0.06 


0.07 


HD107168 


A8m/kA5hA5mF0 


-0.65 


0.08 


-0.40 


0.04 


1.06 


0.03 


0.49 


0.18 


0.77 


0.12 


HD 109307 


A4Vm/A3IV-V 


-0.36 


0.01 


-0.31 


0.12 


0.65 


0.45 


-0.07 


0.29 


0.23 


0.18 


HD 108642 


A2m/kA2hA7mA7 


-0.75 


0.16 


-0.85 


0.12 


0.13 


0.14 


0.25 


0.19 


0.09 


0.09 


HD 107276 


Am/Ka5mA7 


-0.13 


0.27 


-0.08 


0.04 


0.51 


0.08 


0.25 


0.18 


0.21 


0.18 


HD108486 


AmkA3hA5mA7 


-0.69 


0.20 


-0.86 


0.04 


0.60 


0.34 


-0.02 


0.18 


0.02 


0.21 


HD 106999 


Am 


-0.15 


0.10 


-0.02 


0.04 


1.04 


0.11 


0.40 


0.18 


0.26 


0.09 


HD107513 


Am/kA7hF0mF0 


-0.36 


0.10 


-0.15 


0.12 


0.10 


0.10 


0.22 


0.16 


0.12 


0.12 


HD 


SpT 


Call 


a-ca 


Sell 




Till 


CTjj 


VII 


cry 


CrII 


(Tcr 



HD107966 


A3V/A3IV 


-0.08 


0.08 


-0.13 


0.08 


-0.10 


0.05 






0.08 


0.09 


HD108382 


A4V/A3IV 


-0.19 


0.04 


0.12 


0.26 


-0.09 


0.17 






0.24 


0.25 


HD 106887 


A4m 


-0.11 


0.21 


0.03 


0.34 


0.12 


0.15 






0.13 


0.12 


HD107655 


AOV 






-0.20 


0.11 


0.02 


0.08 






0.56 


0.09 


HD107168 


A8m/kA5hA5mF0 


0.23 


0.16 


-0.18 


0.10 


0.48 


0.08 






0.48 


0.20 


HD 109307 


A4Vm/A3IV-V 


0.14 


0.06 


0.10 


0.11 


-0.07 


0.10 






0.03 


0.05 


HD 108642 


A2m/kA2hA7mA7 


-0.33 


0.18 


-1.17 


0.20 


-0.18 


0.07 


0.86 


0.13 


0.18 


0.05 


HD 107276 


Am/Ka5mA7 


-0.24 


0.16 


-0.06 


0.23 


0.04 


0.22 


0.90 


0.10 


-0.08 


0.11 


HD 108486 


AmkA3hA5mA7 


-0.22 


0.15 


-0.50 


0.31 


-0.20 


0.13 


0.70 


0.29 


0.08 


0.15 


HD 106999 


Am 


0.10 


0.16 


0.41 


0.27 


0.14 


0.21 


0.69 


0.29 


0.03 


0.20 


HD107513 


Am/kA7hF0mF0 






-0.25 


0.10 


0.01 


0.07 


0.65 


0.32 


0.09 


0.09 


HD 


SpT 


Mnl 


O" Mil 


Fell 


CTpe 


Co 




Nil 


CTNi 


SrII 


a-sr 



HD107966 


A3V/A3IV 


-0.35 


0.18 


-0.13 


0.05 


0.65 


0.38 


-0.18 


0.07 


-0.28 


0.33 


HD108382 


A4V/A3IV 


-0.55 


0.12 


-0.14 


0.10 






-0.14 


0.16 


-0.26 


0.02 


HD106887 


A4m 


0.02 


0.21 


0.21 


0.15 






0.37 


0.11 


0.59 


0.08 


HD107655 


AOV 






0.08 


0.05 






0.77 


0.23 


-0.23 


0.20 


HD107168 


A8m/kA5hA5mF0 


0.19 


0.12 


0.39 


0.10 






0.60 


0.10 


0.79 


0.31 


HD 109307 


A4Vm/A3IV-V 


-0.04 


0.17 


0.05 


0.09 


-0.06 


0.35 


0.12 


0.08 


0.46 


0.10 


HD 108642 


A2m/kA2hA7mA7 


-0.04 


0.10 


0.16 


0.11 


0.29 


0.35 


0.41 


0.08 


0.50 


0.13 


HD 107276 


Am/Ka5mA7 






0.03 


0.23 






-0.14 


0.10 


-0.22 


0.19 


HD108486 


AmkA3hA5mA7 






0.20 


0.09 






0.21 


0.21 


0.76 


0.02 


HD 106999 


Am 






0.08 


0.13 






0.23 


0.20 


0.47 


0.07 


HD107513 


Am/kA7hF0mF0 






-0.02 


0.07 






-0.23 


0.10 


-0.06 


0.19 


HD 


SpT 


YII 


ay 


Zrll 




Ball 













HD107966 


A3V/A3IV 


0.00 


0.13 


0.38 


0.12 


0.04 


0.26 


HD108382 


A4V/A3IV 


0.25 


0.21 


0.16 


0.17 


-0.24 


0.03 


HD106887 


A4m 


0.78 


0.16 


0.61 


0.15 


1.40 


0.26 


HD107655 


AOV 


0.80 


0.15 


0.78 


0.12 


0.72 


0.32 


HD107168 


A8m/kA5hA5mF0 


0.97 


0.12 


0.93 


0.13 






HD 109307 


A4Vm/A3IV-V 


0.56 


0.07 


0.60 


0.11 


1.17 


0.23 


HD 108642 


A2m/kA2hA7mA7 


0.88 


0.12 


0.75 


0.12 


1.79 


0.24 


HD 107276 


Am/Ka5mA7 


0.12 


0.10 


-0.15 


0.20 


0.44 


0.21 


HD108486 


AmkA3hA5mA7 


0.56 


0.12 


0.67 


0.13 


1.57 


0.21 


HD 106999 


Am 


0.27 


0.11 


0.49 


0.08 


0.78 


0.07 


HD107513 


Am/kA7hF0mF0 


0.08 


0.11 


0.13 


0.20 


0.53 


0.15 



This can be established independently of errors in the absolute 
values of the oscillator strengths, since all stars will be affected 
in the same manner 



4.2.1 . The light elements 
Behaviour of carbon: 

Seven lines of quality B of C I have been synthesized. In a graph 
of [-^l versus Teff, this element displays a different behaviour in 
F stars and A stars. The F stars show very little scatter around 
their mean iron abundance (about -0.01 dex) whereas the A and 
Am stars exhibit a pronounced spread in abundances of about 
0.75 dex which is much larger than the estimated typical uncer- 
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Table 5. Abundances relative to hydrogen and to the solar value, [^] = log(^)* - log(^)o for the F stars. The HD numbers in 
italics are those for which the uncertainties have been calculated as explained in Appendix A. For the others, the quantities labeled 
as cr are standard deviations. 



HD SpT CI o-c OI o-Q Nal 0-^0 Mgll(Mgl) cr„g Sill aSi 



rlUlUoWJ 


r3 V 


-U.U4 


n no 
O.OO 


-U.il 


n 1 c 
U.lj 


-U.U4 


0.04 


0.30(0.02) 


ATI /n nn\ 
U./1(U.U9) 


A 1 A 


A 1 A 


HD106293 


F5V 


-0.07 


0.27 






0.01 


0.18 


0.45(0.10) 


0.18(0.06) 


0.26 


0.11 


HD 106691 


F5IV 


-0.07 


0.16 


-0.26 


0.15 


-0.06 


0.20 


0.30(0.10) 


0.18(0.07) 


0.12 


0.08 


HD 106946 


F2V 


-0.02 


0.20 


0.02 


0.15 


0.09 


0.14 


0.45(0.13) 


0.18(0.04) 


0.26 


0.07 


HD107611 


F6V 


0.05 


0.10 






-0.10 


0.04 


0.23(0.10) 


0.18(0.15) 


0.22 


0.12 


HD 1095 30 


F2V 


0.08 


0.09 


0.10 


0.14 


0.16 


0.12 


0.39(-0.01) 


0.11(0.13) 


0.13 


0.28 


HD107877 


F6 










-0.13 


0.03 


0.40 


0.18 






HD108154 


F5 










-0.13 


0.21 


0.15 


0.18 






HD 108226 


F5 










0.04 


0.23 


0.22 


0.18 






HD108976 


F6 V 










0.02 


0.20 


0.13 


0.18 






HD 109069 


FO V 














0.41 


0.18 






HD 


SpT 


Call 


0"Cfl 


Sell 


0"5c 


Till 


CTri 


Vll 


CTv 


Crll 


(Tcr 
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tainty of about 0. 15 dex on [^]. Also all A stars display deficien- 
cies in carbon. There seems thus to be real star-to-star variation 
in [-^l- Another way to quantify the relative dispersion of A and 
Am stars with respect to the F stars is to compute the mean abun- 
dance for carbon and the associated dispersion for each group of 
stars, (Ta and (Tp. For A and F stars, these mean abundances and 
dispersions for carbon and all other chemical elements are col- 
lected in Table |6] For carbon, the dispersion cta for the A stars 
is about 4 times higher than for F stars. 

The behaviour of carbon with respect to iron is not clear when 



[■^] is displayed versus [^]. As previously found by Hill (1995), 
an anticorrelation appears more clearly when [-^l is displayed 

against {^~\ for Am and normal A stars (Figure |2ll. We find a 
slope of - 1 .74 + 0.48 , not very different from the value obtained 
by Hill (1995) (-1.24 + 0.31 for 15 A stai's). 
We have checked whether the C lines studied here might 
be affected by non-LTE effects especially for A stars. Non- 
LTE abundance corrections of carbon have been calculated by 
|Rentzsch-Holm (1996)| for a set of main sequence stars rang- 
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Fig. 6. Abundance patterns for the "normal" A (a). Am (b) and F stars of Coma Berenices cluster. A maximum +0.30 dex error bar 
is displayed. The horizontal dashed line represent the solar composition. 



ing in effective temperatures from 7000 K to 12000 K, sur- 
face gravities from log^ = 3.5 to 4.5 dex and metallicities 
from [^] = -0.5 to 1 dex. For effective temperatures below 
10000 K, the non-LTE abundance corrections were found to be 
always negative. Only three lines of our list were studied by 
|Rentzsch-Holm (1996)t ^5052.17, /15380.34 and ^6587.61 A. 
According to her Figures 7 and 8, the corrections for these lines, 
whose equivalent width are always less than 100 mA, should be 
in the range -0.10 to -0.25 dex for the temperatures and metal- 
licities of the stars in our sample, making these stars even more 
deficient in carbon. These corrections do not appear to diminish 
the star-to-star scatter in [C/H]. 



Oxygen 

The derived oxygen abundances are the weighted means of 16 
quality B and 6 quality C+ lines. In a graph of [^] versus Teff 
(Figure |7]i, the F stars abundances are slightly scattered around 
their mean value. All A and Am stars exhibit underabundances 
of oxygen and their scatter is again larger than that of the F stars. 
The total spread in oxygen abundance for the A and Am stars is 
about 0.8 dex, significantly larger than the maximum uncertainty 
of 0.12 dex for A stars. This again suggests real star-to-star vari- 
ations in [O/H] . 

In a graph of [^] versus [jf], oxygen seems to be only very 
loosely anticorrelated with iron. The anticorrelation appears 
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Fig. 7. Left panel: Abundance of carbon, oxygen and sodium versus effective temperature. The dotted line corresponds to the solar 
value and the dashed one to the mean value determined for the F stars of the cluster. Right panel: [C/Fe], [O/Fe] and [Na/H] versus 
[Fe/H]. The filled dots correspond to normal A stars, the filled squares correspond to Am stars and the filled diamonds correspond 
to F stars. In the plot representing [Na/H] versus [Fe/H], the dashed line corresponds to the solar [Na/Fe] ratio. The error bars in the 
right panel represent the mean uncertainties for the displayed abundances. 



more clearly in [■!;] v ersus [■^] with a slope of -2.30 + 0.53. 
Przybilla et al. (2000)| calculated the line formation of most of 
the O I lines analysed here in A star atmospheres using the codes 
DETAIL and SURFACE (Giddings 1981). For stars with effec- 
tive temperatures less than 10000 K and \ogg around 4.0, these 
corrections always remain less than -0.03 dex, well below the es- 
timated uncertainty. The observed anticorrelation of [O/Fe] ver- 
sus [Fe/H] should therefore remain after non-LTE corrections 
have been made. The large star-to-star variations in [O/H] would 
not be affected by these corrections either. 



Sodium 

Sodium abundances were derived from 9 lines of quality A to C. 
The sodium abundances for the F stars display very little scat- 
ter around their mean value (-0.01 dex), very close to solar All 
A and Am stars display overabundances of Na. For these stars, 
the spread in [^] is about 0.9 dex, much larger than the typical 
uncertainty (0.10 dex). This suggests there are real star-to-star 
variations in [^]. The sodium abundance does not appear to be 
correlated to the iron abundance: in Figure |2l about half of the 
data lie above the line of solar ratio while the other half are scat- 
tered around it. 

Quantitative information on non-LTE corrections for Na abun- 
dances in A stars is scarce. Non-LTE corrections have been 
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Fig. 8. Left panel: Abundance of magnesium, silicon and calcium versus effective temperature. The dotted line represents the solar 
value and the dashed one represents the mean abundance of F stars. Right panel: [Mg/H], [Si/Fe] and [Ca/H] versus [Fe/H]. The 
symbols are the same as in Figure]?] The dashed lines represent the solar ratios. 



calculated by ]Bikmaev et al. (2002)] for two cool A8IV-V stars 
whose fundamental parameters are close to the coolest A star of 
our sample. The largest corrections amount to -0.35 dex and oc- 
cur for the two Unes /15889.95 A and 5895.92 A. In the absence 
of calculations for hotter stars with various [■^], it is difficult to 
predict what the non-LTE corrections would be for the A and 
Am stars studied here. Provided the corrections would be small 
for F stars, lowering the A and Am star abundances by about 
0.35 dex would certainly improve the correlation of [^] with 

Magnesium 

Seven lines of Mg II of quality B to D were used for the A stars 
and 5 lines of Mg I of qualities B to C for the F stars. In graphs 



of [-^] versus Tgff (Figure ]8j, the F stars and the A and Am 
stars display fairly large and comparable scatter The maximum 
spread in [■^] is about 0.30 dex and 0.55 dex for the F and the 
A and Am stars respectively. The maximum uncertainty in [■^] 
being about 0.18 dex, there does not appear to be significant star- 
to-star variations. Almost all stars exhibit LTE overabundances. 
We noticed that the Mg II A 448 1 A triplet systematicaly yields 
higher abundances than other Mg II lines. Excluding this triplet 
reduces the abundances by about -0.30 dex for the A and F stars. 
The [^] abundances (derived from all lines excluding the Mg 
II triplet) do not appear to be correlated to [■^]. Most data lie 
above the line of the solar ratio [■^]. In their analysis. Hill & 
Landstreet (1993) found that the ratio [j^] runs parallel to and 
above the line of solar [y^]. They attribute this to the use of the 
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Fig. 9. Left panel: Abundance of scandium, titanium and vanadium versus effective temperature. The dotted line represents the solar 
value and the dashed one represents the mean abundance of F stars. Right panel: [Sc/H], [Ti/Fe] and [V/H] versus [Fe/H]. The 
symbols are the same as in Figure]?] The dashed lines represent the solar ratios. 



very strong Mg II A 4481 A triplet (the only line they analysed). 
Przybilla et al. (20 01)]computed a model atom for non-LTE line 
formation for neutral and singly-ionized magnesium to evaluate 
non-LTE corrections. These corrections turn out to be small for 
Mgll except for the features at AA 4481 and 7877-96 A. For Vega 
( AOV), they found a correction of -0.2 1 dex for the 448 1 A triplet 
which confirms our finding that this line systematically yields 
higher abundances. 

Silicon 

Nineteen lines of Si II of quality C to E have been synthesized. 
For the F stars, the silicon abundances hardly show any scatter 
around their mean value, -1-0.18 dex, slightly above solar Almost 
all A stars are also overabundant in silicon and deviate very little 



from this mean value except for the Am star HD 107168 (Figure 
]8]l. Thus there does not seem to be significant star-to-star varia- 
tion in [^]. 

The silicon abundance does not show any convincing correlation 
with the iron abundance. Most of our data fall slightly above the 
line representing a solar silicon to iron ratio. Note that Hill & 
Landstreet (1993) did find a tight correlation between [^] and 
[■^]. All their data fall close to the line representing a solar sil- 
icon to iron ratio. The lines we analysed differ from theirs and 
lead to overabundances possibly because of incorrect oscillator 
strengths or non-LTE effects. 
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Fig. 10. Left panel: Abundance of chromium, manganese and iron versus effective temperature. The dotted line represents the solar 
value and the dashed one represents the mean abundance of F stars. Right panel: [Cr/H] and [Mn/H] versus [Fe/H]. The symbols 
are the same as in Figure]?] The dashed lines represent the solar ratios. 



Calcium 

Twelve lines of Ca II of quality C and D were used to derive 
the calcium abundance for all A and Am stars and for only two 
F stars. The abundance of the F stars is subsolar (< [^] >~ 
-0.2 dex). The two normal A stars are slightly deficient in cal- 
cium and the Am stars show modest deviations (both over and 
underabundances) around the solar abundance. Their maximum 
spread in [Ca/H] is 0.50 dex which is only marginally signif- 
icant compared to the maximum estimated uncertainty (about 
0.20 dex). 

In a graph displaying the calcium abundance versus that of iron 
(FigurelSJ, the Am stars, as expected, can easily be discriminated 
from the normal A and F stars. They all fall to the right in a re- 
gion characterized by overabundances of iron and they are either 
Ca-deficient or Ca-rich. For the normal A and F stars, [^] has 



a fairly uniform value (-0.20) dex and does not appear to vary 
with 1^1 

Scandium 

Eleven lines of quaUty D of Sc II were used to derive the scan- 
dium abundance. The scandium abundances of F stars are scat- 
tered very little around the solar value. About half of the A and 
Am stars are deficient in scandium (3 Am stars are close to solar 
while 4 exhibit deficiencies ranging from -0.2 dex to -1.2 dex). 
The total spread in [^] for these stars is about 1.20 dex, signifi- 
cantly larger than the typical uncertainty of 0. 10 dex. There thus 
seems to be real star-to-star variation in [^]. 
As for calcium, scandium does not exhibit any clear correlation 
or anticorrelation with respect to iron in the diagram of [^] ver- 



16 



M. Gebran, R. Monier and O. Richard: Chemical composition of A and F dwarf members of the Coma Berenices open cluster 



-1 



-2 



-2 





# Normal A 'stars 


: [co/H] 


■ Am stars 
4 F stars 


♦ T 

: ♦ 






Solar value 

— — Mean value (F stars) 


1 1 1 1 1 1 1 1 

: [Ni/H] 


1 1 1 1 1 1 1 1 1 


1 

- 


i 

■ ■ 

! f • 

- 


1 1 1 1 1 1 1 1 

: [Sr/H] 


1 1 1 1 1 1 1 1 1 




T 


L»,:^.^....r...... 


■ • • 







6500 7000 7500 8000 8500 9000 



• ♦ 



0.5 - 







♦ ♦ 



-0.5 



-1 



[Co/H] 

Solar ratio 




Regression for A stars 




ression for A stars _ 



-1.5 



-0.4 



-0,2 



0,2 
[Fe/H] 



0,4 



0,6 



Fig. 11. Left panel: Abundance of cobalt, nickel and strontium versus effective temperature. The dotted line represents the solar 
value and the dashed one represents the mean abundance of F stars. Right panel: [Co/H], [Ni/Fe] and [Sr/H] versus [Fe/H]. The 
symbols are the same as in Figure]?] The dashed lines represent the solar ratios. 



sus [^] (Figure ]9]l. The scandium abundances of the F and nor- 
mal A stars, which are all very close to solar, do not depend on 
[■^]. Most Am stars lie in the lower right part of the diagram: 
they are all iron-rich and 4 out of 7 are deficient in scandium. 

4.2.2. The iron peak elements 

Titanium, vanadium, cliromium and manganese 

Twenty six lines of Ti 11, most of them of quality D, have been 
synthesized. The oscillator strengths for this element are not as 
secure as for iron and chromium. One should not expect these 
abundances to be too reliable on an absolute scale. The F stars 
show little scatter around their mean [jj-] value, about -1-0.12 
dex, slightly above solar. The A and Am stars are more scat- 



tered around this mean value. The spread in [^] is about 0.68 
dex which is larger than the typical uncertainty of about 0.08 
dex. There seems to be real star-to-star variation for this element 
(Figure]9]l. The titanium abundance appears to be loosely corre- 
lated to that of iron (correlation coefficient 0.66) for normal F 
and A stars. The [j^] ratios of these normal stars appear to be 
close to or slightly higher than solar as found by Lemke (1989) 
for a sample of 16 normal A stars. Hill & Landstreet (1993) 
found a strong correlation of [Ti/H] with [Fe/H] (correlation co- 
efficient 0.95), their [j^] ratios being slightly above the solar 
value. 

Nine lines of V 11 were synthesized whose uncertainties are not 
specified. These lines are intrinsically weak in all spectra. The 
abundances derived for this element must be taken with cau- 
tion. The vanadium abundances for the F stars are fairly scat- 
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Fig. 12. Left panel: Abundance of yttrium, zirconium and barium versus effective temperature. The dotted line represents the solar 
value and the dashed one represents the mean abundance of F stars. Right panel: [Y/H], [Zr/Fe] and [Ba/H] versus [Fe/H]. The 
symbols are the same as in Figure|2l The dashed lines represent the solar ratios. 



tered around their mean value, +0.50 dex. Vanadium was not 
detected in the two normal A stars and is found to have large 
overabundances for 4 Am stars. There is no convincing correla- 
tion between [^] and [^] (Figure |9]l. Most of the data fot the 
normal A stars and the F stars fall above the line representing 
the solar [j^] ratio which may be due to poorly determined os- 
cillator strengths. Hill & Landstreet (1993) found a correlation 
of [■^] with [^] using lines other than ours; their ratios [j^] are 
also significantly above solar. 

Eleven lines of Cr II of quality D have been synthesized. As for 
titanium, the A stars are only a little more scattered than the F 
stars. The total spread in [^] for the A and Am stars is about 
0.56 dex (0.36 dex for the F stars) while the maximum uncer- 
tainty on [^] is about 0.20. Should the data for HD 107 168 be re- 
moved, the spread in [£^] drops to 0.30 dex and is not significant. 



The evidence for star-to-star variations in [^] is therefore rather 
weak. The chromium abundance appears to be loosely correlated 
with that of iron (correlation factor = 0.85). Most of the data fot 
the normal A stars and the F stars fall above the line represent- 
ing the solar [^] ratio (Figure[TOll. Hill & Landstreet (1993) also 
found a correlation of [^] with [■^], their ratios [^] being only 
marginally above solar. 

Twenty lines of Mn I of quality B to C+ were synthesized. The F 
stars show little scatter around the solar value. The total spread 
for the A and Am stars is about 0.6 dex which is larger than 
the maximum uncertainty of 0.2 dex, suggesting real star-to-star 
variation. The two normal A stars are deficient in manganese. 
The manganese abundance appears to be well correlated with 
that of iron (coiTelation factor = 0.94). Most of the ratios [^] 
run parallel to and below the line representing the solar ratio 
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(Figure[TOll. Hill & Landstreet (1993) found a correlation of [^] 
with [jf-] using lines other than ours; their ratios [yj] are only 
marginally above solar. 



Iron, cobalt and nickel: 

Twenty seven lines of Fe 11 of quality C to E were synthesized. 
In a graph [^] versus T^ff (Figure [TO] i. the F stars are only 
slightly scattered around their mean value, +0.07 dex, which 
is slightly higher than Friel & Boesgaard (1992) These authors 
found < [^] >- -0.05 + 0.03 dex, based on the analysis of 



Table 6. Mean abundances and their respective dispersion for 
Fe, C, O, Ni, Sc, Si, Ba, Y, Zr and Sr for the F, normal A and Am 
stars. 



the equivalent widths of a few Fe 1 lines for 14 F stars of the 
cluster Our usage of different techniques (model atmospheres 
and line synthesis) and different lines (Fe 11) probably accounts 
for the difference in iron abundance. The total spread in [^] for 
the A and Am stars is about 0.53 dex which is larger than the 
maximum estimated uncertainty (0.10 dex). There are thus real 
star-to-star variations in [^]- 

For cobalt, our analysis is based on 12 lines of Co I, whose er- 
rors in the oscillator strengths are unknown. Most of these lines 
are weak and often are blended with lines whose atomic param- 
eters are not necessarily accurately known. We therefore do not 
expect these abundances to be reliable. In a graph of [^] versus 
Teff (Figure [TT]), the cobalt abundances are much more scattered 
for the F stars than for any other chemical element; this scatter is 
probably largely due to blending species. There is only one data 
point for the normal A stars. We feel that the abundances are not 
reliable enough to claim to find real star-to-star variations in [^] 

and not to find a correlation with [^]- 

Fifty six lines of Ni 1 of quality C+ to D have been synthesized. 
In the graph [jj-] versus Tgff, nickel behaves in a similar manner 
as iron (Figure fTTTi. The F stars are fairly well grouped around 
their mean abundance {+ 0.09 dex), their maximum spread be- 
ing 0.25 dex. The A and Am stars scatter over 0.83 dex whereas 
their maximum uncertainty is about 0.20 dex. There seems to 
be real star-to-star variations for [jj-]. The nickel abundance is 
tightly correlated with that of iron (correlation coefficient 0.92). 
Most of the ratios [|^] are close to solar 



4.2.3. The heavy elements 

Strontium, yttrium, zirconium and barium: 

Our study of strontium is based on only two lines at 4077.71 
A and 4215.52 A whose accuracies are not specified in NIST. 
Errors in gf values may result in a considerable zero point shift 
with respect to the Sun. We therefore do not attach too much 
significance to apparently large absolute over-or underabun- 
dances. The F star strontium abundances show little scatter 
around their mean value, + 0.15 dex. The A and Am stars 
display a much larger spread in abundances, about 1.07 dex, 
significantly larger than the maximum uncertainty (0.30 dex) 
(Figure [TTli. The strontium abundance is correlated quite closely 
to the iron abundance (correlation coefficient 0.90). Strontium 
abundances vary more rapidly than those of iron as found by 
Lemke(1990). 

Five lines of yttrium whose accuracies are unknown were 
synthesized. The same holds for zirconium. Yttrium and zirco- 
nium are found to be nearly solar in F stars (Figure [T2]l. These 
elements are overabundant by about 0.50 dex for 6 of the Am 
stars and by about 0.40 dex in the normal A stars. The A and 
Am stars display fairly large spreads in abundances of Y and 
Zr, about 1.0 dex, much larger than the associated uncertainties. 



Elements 


F stars crp 


A stars cta 


[^] 


-0.01 0.06 


-0.48 0.24 




-0.11 0.18 


-0.34 0.31 




0.18 0.07 


0.22 0.20 


r— 1 
— — 


-0.02 0.04 


-0.23 0.38 


r— 1 


0.07 0.09 


-0.14 0.16 


— — 


0.09 0.11 


0.11 0.28 




0.15 0.08 


0.25 0.43 




0.07 0.08 


0.47 0.34 




0.04 0.34 


0.46 0.33 




0.58 0.17 


0.84 0.70 



indicating real star-to-star variations in [-^l and [jf]- 
Five lines of quality B of barium have been synthesized. This 
element is found to be overabundant in all A, Am and F stars 
by large amounts (up to 1.10 dex). The spread for [^] is about 
1.8 dex, much larger than the typical uncertainties, indicating 
real star-to-star variations in [^] (Figure [T2]i. However, these 
overabundances may reflect a non-LTE effect, namely an 
overionization of barium in A stars. In Vega, the non-LTE Ba 
II abundances are lower than the LTE abundances by about 
0.30 dex as demonstrated by |Gigas (1988)] and [Lemke (1990)| 
using the KIEL code. Non-LTE corrections for Ba should be 
sensitive to [Fe/H] which can differ much from star to star. 
Detailed calculations for the respective fundamental parameters 
and iron abundances of each star should therefore be performed . 



4.3. Microturbulent velocities 

A byproduct of this analysis has been to determine microturbu- 
lent velocities for each star Figure [13] displays the derived ^, 
versus effective temperature. The overall variation agrees well 
with that found by Coupry & Burkhart (1992) who found that ^, 
varies from km-s"' for late B-type stars, up to about 3 km-s"' 
for mid-A type stars down to around 2 km-s ' for early F-type 
stars. Gray (2001) also found that ^, varies diminishes 3 km-s"' 
for mid-A type stars to about 1 km-s"' for solar-type stars. 



4.4. Rotationai velocities 

The inferred rotational velicities VfSin/ are lower than 102 
km-s"'. None of the found abundances appears to depend on 
Ve sin /. The profiles of [X/H] versus Vf sin / are flat as shown 
for instance for iron in Figure [14] This result is not surprising. 
Detailed calculations of diffusion in the presence of meridional 
circulation carried out by Charbonneau & Michaud (1991)| re- 
vealed that, in stars rotating at less than VcSin;= 100 km-s~', 
meridional circulation has little influence on chemical separation 
once the helium superficial convective zone has disappeared. 
Accordingly, the abundances should not present any positive nor 
negative trend (slope) with v^. sin / in this velocity regime. The 
absence of fast rotators in Coma Berenices (ie. Vf,sin/ > 100 
km-s"') prevents us from investigating the behaviour of the var- 
ious abundances with rotational velocity above that velocity. 
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(km/s) 



Normal A stars 
Am stars 
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Fig. 13. Variation of the derived microturbulence velocities with 
effective temperature. 



[Fe/H] 



Normal A stars 
Am sl.ars 



Fig. 14. Abundance of iron versus Ve sin / for A,Am and F stars. 



5. Astrophysical implications 

The most important result of this study is the evidence of large 
star-to-star abundance variations for A stars in Coma Berenices. 
These stars appear to display much larger star-to-star variations 
in their abundances than the F stars do for the following chemical 
elements: C, O, Na, Sc, Ti, Mn, Fe, Ni, Sr, Y, Zr and Ba. In con- 
trast, the abundances of Mg, Si, Ca and Cr do not show signif- 
icant star-to-star variations. For the Hyades, Varenne & Monier 
(1999) found a similar behaviour for the abundances of C, O, 
Na, Sc, Fe and Ni. Monier (2005) also found large star-to-star 
variations in [Fe/H], [Ni/H] and [Si/H] for several A stars of the 
Uma group. 

We theorize that this peculiar behaviour is a signature of the 
occurence of transport processes competing with radiative dif- 
fusion (eg. rotational mixing in radiative zone, Zahn 2005). 
Indeed, if radiative diffusion was the only process at work, 
the microscopic diffusion velocity of a given chemical element 



should be the same in stars of similar effective temperatures 
and surface gravities. We would therefore expect similar surface 
compositions for stars of similar fundamental parameters. 

5.1. Self consistent modelling 

We have compared the derived abundances to the predictions of 
recent evolutionary models, calculated with the Montreal code 
using slightly different assumptions and physics. This code treats 
radiative diffusion in detail and allows the inclusion of turbulent 
diffusion. Schatzman (1969) proposed two possible physical ori- 
gins of turbulent diffusion: loss of angular momentum while the 
star descends towards the Main Sequence or meridional circula- 
tion. 



5.1.1. The F-type stars 

For the F stars, the found abundances have been compared to the 
predictions of Turcotte et al's (1998) evolutionary models at the 
age of Coma Berenices. Their models predict the evolution of 
the abundances for an F star of a given mass consistently with the 
internal structure. They include the effects of gravitational sed- 
imentation and of radiative diffusion for 28 chemical elements 
(Z < 28) but no macroscopic mixing (wind, accretion, merid- 
ional circulation, turbulence). These models are relevant for F 
stars with masses in the range I.IMq to l.SM© from the pre- 
main sequence state up to hydrogen core exhaustion. The impact 
of the abundance variations with time on the structure of the star 
is taken into account. Monochromatic OPAL opacities for each 
element are used to recalculate the opacity corresponding to the 
abundances and local conditions during the evolution. Atomic 
diffusion has an important effect on the opacities for stars more 
massive than 1.3Mq. In these objects, the abundances of Fe and 
other iron-peak elements were found to substantially vary with 
time. 

For the light elements C and O, the predicted large underabun- 
dances by Turcotte et al. (1998) for stars with T^ff > 6500 K 
are not observed in our data. For these elements we find that the 
observed abundances in F stars are rather constant with a solar 
value for C (< [§] >= -0.01+0.06 dex) and marginally under- 
abundant for O (< [§] >= -0.11+0.18 dex). In contrast, the 
predicted solar Na abundances for stars with 5950 < Tejj < 
6700 K match well the observed solar abundances of F stars. For 
the F stars of I.AMq, the predicted Mg and Si underabundances 
are not seen either in our data. For the iron peak elements Fe and 
Ni, the predicted overabundances due to microscopic diffusion 
are not observed either. 



5.1.2. The A stars 

For A stars, we have compared the found abundances with 
the predictions of the recent models of [Richer et al. (2000)| In 
these new models, the effect of atomic and turbulent diffusion 
were calculated for stars of 1.45-3 Mq. They showed that the 
superficial abundances of the 28 species calculated in their 
models depend, in a star of a given mass, on essentially the 
initial metallicity and the deph of the zone mixed by turbulence. 
Figures 10 and 11 from |Richer et al. (2000)] represent the 
variations of surface abundances with time of the 28 elements 
for stars of 2-2.5-3 Mq. For the 2 Mq model, [Richer et al. (2000)] 
found that the ratio of the abundances of carbon and oxygen 
at 450 Myr to their initial values decrease with time. We 
observed the same trend for C and O in Am stars using the mean 
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abundances of C and O in the F stars as initial abundances of 
these elements in the cluster. 

In addition, we have specifically computed a series of mod- 
els for two Am stars of similar effective temperatures but dif- 
ferent rotational velocities: HD108642, a slow rotator (V(.sin/ 
= 9 km-s"') and HD106887, a faster rotator (v^sin/ = 82 
km-s"'). Masses for these two stars should be in the range 
1.8 to 2.0 Mq. Seven models for a 1.8 Mq mass star hav- 
ing different turbulent coefficients, Dj (equation 1 in Richer 
et al. 2000) at the age of the Coma were computed. The 
adopted initial and homogeneous abundances for these mod- 
els are collected in Table [T] The models use an Eggleton- 
Faulkner-Flannery equation of state ( Eggleton et al. 1973) in- 
cluding Coulomb correction on the pressure (labeled as CEFF 
models) (see also [Christensen-Dalsgaard & Daeppen 1992) l. 
The nuclear energy generation follows the prescriptions of 
[Bahcall & Pinsonneault (1992)| These models take into account 
gravitational settling, thermal diffusion and radiative accelera- 
tions. The detailed treatment of atomic diffusion is described 
in|Turcotte et al. (1998) and the radiative accelerations are from 



Turcotte et al. (1998) with correction for redistribution from 
Gonzalez etal. (1995), and [LeBlanc et al. (2000)| These mod- 



els are self consistent as the Rosseland opacity and radia- 
tive accelerations are recomputed at each time step in each 
layer for the exact local chemical composition using OPAL 
monochromatic opacities for 24 elements. Convection and semi- 
convection are modeled as diffusion processes as described in 
[Richer et al. (2000)| 



and Richard etal. (2001) 



^ , The initial metal- 

licity is taken from 
The models are compared to the abundance patterns of 
HD108642 in Figure [B] and for HD106887 in Figure[l6] None 
of the models reproduces entirely the characteristic abundance 
pattern, ie. marked underabundances of light elements and over- 
abundances of iron-peak elements. They basically differ by the 
amount of turbulent diffusion included. The model that best ap- 
proaches the abundance patterns for elements with Z>20 for 
these two stars is the model labeled as 1.8T5.3D200K-3 follow- 
ing the syntax of Table 1 of Richer et al. (2000)] In this model, 
the turbulent diffusion coefficient Dt varies with density as: 



Dj ^ ajD{He)o{f)" 



where n=3, D(He)o is the atomic diffusion coefficient of He at 
the density po - piTo) in the iron convection zone (see equation 
1 of Richer et al. 2000, here w = 200 ■ 10^ and log Tq = 5.3). 
However this model predicts abundances that are almost 0.8 dex 
too large for C and O but it comes close to the abundances of 
Na, Mg and Si and follows the overabundances of iron-peak 
and heavier elements. Conversely, models with less turbulent 
diffusion 1.8T5.3D25K-3 and 1.8T5.3D500-3 roughly account 
for the abundances of elements with Z<15 but predict too large 
overabundances of iron-peak elements. 

Part of the discrepancy between observed and theoretical abun- 
dance patterns could be due to non-LTE effects. Correcting the 
magnesium abundance for the magnesium triplet at /14481 A 
reduces it by -0.2 dex and provides a better agreement with 
the T5.3D200K-3 model for this element. Lowering the sodium 
abundances, which are likely to be affected by non-LTE effects, 
will also improve the agreement with this model. However, the 
inclusion of competing processes in the models such as rota- 
tional mixing in the radiative zone for A stars, internal waves for 
F stars could also help improve the agreement. 



Table 7. Initial chemical composition 

Element Mass fraction 

H 7.03x10-' 

"'He*"' 2.7995x10-' 

'2C<*> 2.935x10-3 

N 9.000x10-'' 

O 8.189x10-3 

Ne 1.675x10-3 

Na 3.396x10-' 

Mg 6.377x10-"* 

Al 5.519x10-5 

Si 6.878x10-" 

P 5.944x10-" 

S 3.592x10-" 

CI 7.642x10-* 

Ar 9.170x10-5 

K 3.396x10-' 

Ca 6.368x10"' 

Ti 3.396x10-* 

Cr 1.698x10-' 

Mn 9.340x10-" 

Fe 1.219x10-3 

Ni 7.557x10-' 

(a) 'He = 5.000 x 10-' 

(b) '3Cis l%of '^C 



1.5 



0.5 



-1.0 



-1.5 



1 .8 Af„ T5.3D200k-3 451 Myr, T„= 8190 K, logg= 4.19 ■ 
1.8MgT5.3D25k-3 451 Myr, T„= 8164 K, logg= 4.19 i.. 
r.8Mj,T5.3D500-3 445 Myr, T„= 8113 K, logg= 4.18 »- 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 
H Li B N Na Al P CI K Mn 

He Be C O Ne Mg Si S Ar Ca Ti Cr Fe m 



A HD1 08642 




Fig. 15. Comparison of the observed abundance pattern of 
HD 108642 (A2m) with the predictions of three models calcu- 
lated for a mass of 1.8 Mq and different amounts of turbulent 
diffusion. Observed abundances are represented as triangles with 
error bars. Note that the models do not predict the surface abun- 
dances of Sc, V and Co. 



6. Conclusion 

High and medium resolution spectra of 11 A and 11 F star 
members of the Coma Berenices open cluster have been synthe- 
sized in order to determine the abundances of C, O, Na, Mg, Si, 
Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Sr, Y, Zr and Ba. In graphs 
representing the abundance [X/H] versus effective temperature, 
the A stars display abundances that are more scattered around 
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445 Myr, r_„= 8113 K, logg= 4.18 
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T" 
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l\iln 
Cr Fe 




Fig. 16. Comparison of the observed abundance pattern of 
HD106887 (A4m) with the same models. 



reproduce the observed abundance patterns and the large scatter 
of abundances of several elements in A stars. 



Appendix A: Determination of uncertainties 

Six major sources are included in the uncertainty determina- 
tions: uncertainty on the effective temperature (crj-|^), on the 
surface gravity (crjog^), on the microturbulent velocity (cr^,), 
on the apparent rotational velocity (cr,,^ sin ,), the oscillator 
strength (crioggy) and the continuum placement (cTconi)- These 
uncertainties are supposed to be independent, so that the total 
uncertainty cr,o,, for a given transition (i) is: 



2 J. 2 J. 2 



(A.l) 



The mean abundance < [^] > is then computed as a weighted 
mean of the individual abundances [X/H], derived for each 
transition (i): 



< [^] 
H 



2/(1 /<) 

and the standard deviation, cr,^ is given by: 



(A.2) 



the mean value than the F stars. Large star to star variations are 
detected for A stars for C, O, Na, Sc, Ti, Mn, Fe, Ni, Sr, Y, Zr 
and Ba which we interpret as evidence of transport processes 
competing with radiative diffusion. 

The chemical pattern found for the A and F dwarfs of 
Coma resembles that found for the Hyades (Varenne & 
Monier 1999) and the UMa group (Monier 2005). The 
mean iron abundance derived for the F stars is found to be 
< [^] >= 0. 07 ± 0.09 dex, slightly hig her than the metallicity 
derived by Friel & Boesgaard (1992) The abundances of 



manganese, nickel, strontium and barium are strongly correlated 
with the iron abundance for A and Am stars. The ratios [C/Fe] 
and [O/Fe] appear to be anticorrelated with [Fe/H]. The ratio 
[Ti/Fe] is solar as found by Lemke (1989). 
The Am stars in Coma Berenices are found to be deficient in 
light elements (C and O), they are not all deficient in calcium 
and scandium but are all overabundant in metallic and heavy 
elements (Fe, Co, Ni, Sr, Y, Zr and Ba). The two normal A 
stars have almost solar abundances. The F stars have solar 
abundances for almost all the elements except for Mg, Si, V and 
Ba. 

The abundance patterns predicted by current state of the 
art evolutionary models following the prescriptions of 
[Richer et al. (2000)| have been compared to the observed 
abundance patterns of two Am stars of the cluster, one a 
slow rotator (HD 108642) and one a moderately fast rotator 
(HD106887). These models were calculated with different 
strengths of the turbulent diffusion coefficient. None of the mod- 
els reproduces entirely the characteristic abundance pattern, ie. 
marked underabundances of light elements and overabundances 
of iron-peak and heavier elements. Part of the discrepancy may 
arise from non-LTE effects. In this respect, non-LTE abundance 
determinations for many of the elements analysed here (when 
feasible, ie. when atomic data and model atoms are available) 
are highly desirable. However it is likely that the inclusion of 
competing processes such as rotational mixing in the radiative 
zones (which will vary from star to star) should also help 



1 ^ 

" sd i=l 



crL) 



(A3) 



where N is the number of lines per element. This procedure has 
been applied to 4 A stars (HD107966, HD107655, HD107168 
and HD107513) and 3 F stars (HD106103, HD107611 and 
HD109530). In the fastest rotators, the pseudo-continuum, 
f^' , , we see that regions free of lines in the observed spec- 
tra are actually affected by the rotational broadening of neigh- 
bouring lines. The level of the actual continuum, f^„„,, can be 
recovered by carefully synthesizing the pseudo-continuum win- 
dows, yielding the theoretical flux, F"^ , , and the actual cor- 

responding continuum flux, /^^onr' appropriate velocity 

Ve sin / + A(ve sin /). The observed intensity level in these pseudo- 
continuum windows was then multiplied by {FamJ^'pseudc) '■° 
cover the observed continuum level: 



J CO. 



cont rA 

J pseudo' 



pseudo 



(A.4) 



The maximum and minimum allowed rotational velocity 
V(,sin/ ± A(veSin/) yields 2 ratios (Vg sin imat) and 

pseudo 

pi 

^/°"' (Vf sin /„„■„). The corresponding observed normalized line 

pseudo 

profiles were then used to derive the corresponding changes in 
abundances due to the different locations of the continuum. This 
test was performed in several spectral regions (excluding over- 
laping regions of 2 orders) and yields a maximum uncertainty of 
about 0.07 dex on the abundances. 

For the other stars, the final abundances are averages. The 
errors on the elemental abundances are standard deviation as- 
suming a Gaussian distribution of the abundances derived from 
each line 



N 



(A.5) 
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9 Site - xf 

where x is the mean value of the abundance, N the number of 
lines of the element and a the standard deviation. 
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Table 8. Linelist used for abundance determination. The "A" term is used for a gf accurancy lower than 3%, "B" lower than 10%, "C+" lower 
than 18%, "C" lower than 25%, "D+" lower than 40 %, "D" lower than 50 % and "E" higher than 50 %. If no accurancy is available, we used 
the E {i.e. 50%) value in the un certainties calculations. R eferences are FMW for|Fuhr et al. (1988)| FT P for [Pickering et al. (2002)| NBS for 
Miles & Wiese (1969)| K083 for'K ostyk & Orlova( 1983)1 S L90 for Sigut & Landst reet (1990)| BL48 forlBierinann & Lubeck(1948)' Bal81 for 
Biemont et al.(1981)| MC86 for Ma gazzu & Cowley (1986)1 W85 for Ward (1985) and Kurucz (http://kurucz.harvard.edU/LINELlSTS/GFALL/l 
for the gfall.dat linelist. 
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Table 8. continued. 



Element 
(Ionization) 


AiA) 




Accuracy 


references 


Element 
(Ionization) 


a 

AiA) 




Accuracy 


references 


Till 


4411,0720 


-0,670 


D 


FTP 


Fell 


4233,1720 


-2,000 


C 


FMW 


Till 


4417,7140 


-1,190 




PTP 


Fell 


4258,1540 


-3,400 


D 


FMW 


Till 


4443,8010 


-0,720 


D 


FTP 


Fell 


4273,3260 


-3,258 


D 


FMW 


Till 


4468,4920 


-0,600 


D 


FMW 


Fell 


4296,5720 


-3,010 


D 


FMW 


Till 


4488,3250 


-0,510 


D 


PTP 


Fell 


4385,3870 


-2,570 


D 


FMW 


Till 


4501,2700 


-0,770 


D 


PTP 


Fen 


4416,8300 


-2,600 


D 


FMW 


Till 


4549,6210 


-0,470 


D 


PTP 


Fell 


4472,0900 


-1,791 




Kurucz 


Till 


4571,9710 


-0,320 


D 


PTP 


Fell 


4472,6200 


-2,340 




Kurucz 


Till 


4589,9580 


-1,780 


D 


PTP 


Fell 


4472,9290 


-3,430 




Kurucz 


Till 


4629,2740 


-2,240 




Kurucz 


Fell 


4491,4050 


-2,690 


C 


FMW 


Till 


4657,2000 


-2,240 




Kurucz 


Fell 


4508,2880 


-2,210 


C 


FMW 


Till 


4805,0850 


-1,100 


D 


PTP 


Fell 


4515,3390 


-2,490 


C 


FMW 


Till 


5129,1560 


-1,400 


D 


K083 


Fell 


4520,2240 


-2,600 


C 


FMW 


Till 


5188,6870 


-1,050 


D 


PTP 


Fell 


4522,6340 


-2,030 


C 


FMW 


Till 


5336,7860 


-1,590 


A 


PTP 


Fell 


4541,5240 


-3,050 


C 


FMW 












Fell 


4555,8900 


-2,290 




FMW 


VII 


4475.6700 


-1.440 




Kurucz 


Fell 


4576,3400 


-3,040 




FMW 


VII 


4528.5000 


-0.960 




Kurucz 


Fell 


4582,8350 


-3,100 


C 


FMW 


VII 


4532.1700 


-0.760 




Kurucz 


Fell 


4620,5210 


-3,280 


D 


FMW 


vn 


4538.6200 


-1.800 




Kurucz 


Fell 


4635,3160 


-1,650 


D 


FMW 


VII 


4558.4500 


-0.930 




Kurucz 


Fell 


4656,9810 


-3,630 


E 


FMW 


VII 


4564.5900 


-1.450 




Kurucz 


Fell 


4666,7580 


-3,330 


D 


FMW 


VII 


4577.1300 


-2.140 




Kurucz 


Fell 


4923,9270 


-1,320 


C 


FMW 


VII 


4590.5000 


-0.780 




Kurucz 


Fell 


5197,5770 


-2,100 


C 


FMW 


VII 


4600.1800 


-1.360 




Kurucz 


Fell 


5276,0020 


-1,940 


C 


FMW 












Fell 


5316,6150 


-1,850 


C 


FMW 


CrII 


4558,6500 


-0,660 


D 


Kurucz 


Fell 


5506,1950 


0,950 


D 


FMW 


CrII 


4588,1990 


-0,643 




Kurucz 












CrII 


4592,0490 


-1,217 


D 


FMW 


Col 


4466.8800 


-0.540 




Kurucz 


CrII 


4616,6290 


-1,291 




SL90 


Col 


4469.5400 


-0.330 




Kurucz 


CrII 


4618,8030 


-1,110 


D 


Kurucz 


Col 


4471.5400 


-0.770 




Kurucz 


CrII 


4634,0700 


-0,990 




Kurucz 


Col 


4530.9500 


0.150 




Kurucz 


CrII 


4812,3370 


-1,995 


D 


Kurucz 


Col 


4533.9800 


-0.500 




Kurucz 


CrII 


5237,3290 


-1,160 


D 


FMW 


Col 


4549.6500 


-0.330 




Kurucz 


CrII 


5308,4400 


-1,810 


D 


FMW 


Col 


4565.5800 


-0.220 




Kurucz 


CrII 


5313,5900 


-1,650 


D 


FMW 


Col 


4581.5900 


-0.150 




Kurucz 


CrII 


5502,0670 


-1,990 


D 


FMW 


Col 


4594.6300 


-0.080 




Kurucz 












Col 


4596.8900 


-0.010 




Kurucz 


Mnl 


4451,5860 


0,278 


B 


FMW 


Col 


4625.7600 


-0.370 




Kurucz 


Mnl 


4453,0120 


-0,490 


C+ 


FMW 


Col 


4629.3600 


-0.190 




Kurucz 


Mnl 


4457,0440 


-0,555 


C-H 


FMW 












Mnl 


4458,2540 


0,042 


C-i- 


FMW 


Nil 


4468.4340 


-1.642 




Kurucz 


Mnl 


4461,0790 


-0,380 


C-i- 


FMW 


Nil 


4470.4720 


-0.310 


D 


Kurucz 


Mnl 


4462,0310 


0,320 


C-i- 


FMW 


Nil 


4480.5610 


-1.491 




Kurucz 


Mnl 


4464,6820 


-0,104 


B 


FMW 


Nil 


4490.0490 


-2.108 




Kurucz 


Mnl 


4470,1440 


-0,444 


B 


FMW 


Nil 


4490.5250 


-2.324 




Kurucz 


Mnl 


4472,8060 


-0,583 


B 


FMW 


Na 


4512.9860 


-1.470 


D 


Kurucz 


Mnl 


4490,0900 


-0,521 


B 


FMW 


Nil 


4519.9790 


-2.880 


D-l- 


FMW 


Mnl 


4498,9020 


-0,343 


B 


FMW 


Nil 


4521.3220 


-0.949 




Kurucz 


Mnl 


4502,2130 


-0,344 


B 


FMW 


Nil 


4523.6940 


-1.305 




Kurucz 


Mnl 


4709,7120 


-0,339 


B 


FMW 


Nil 


4528.5260 


-1.127 




Kurucz 


Mnl 


4739,0870 


-0,490 


B 


FMW 


Nil 


4542.2370 


-1.308 




Kurucz 


Mnl 


4754,0420 


-0,085 


B 


FMW 


Nil 


4546.9200 


-0.271 




Kurucz 


Mnl 


4/Di,jiZU 


-U,l JO 


T3 
D 


rMW 




4jji.zi /U 


-U.ooU 


U 


rMW 


Mnl 


4762,3670 


0,426 


B 


FMW 


Nil 


4559.9210 


-1.737 




Kurucz 


Mnl 


4783,4270 


0,042 


B 


FMW 


Nil 


4572.0410 


-0.536 




Kurucz 


Mnl 


4823,5240 


0,144 


B 


FMW 


Nil 


4588.4110 


-0.745 




Kurucz 


Mnl 


5255,3260 


-0,763 


B 


FMW 


Nil 


4592.5220 


-0.370 




Kurucz 












Nil 


4596.3830 


-0.704 




Kurucz 












Nil 


4600.3550 


-0.610 




FMW 
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Table 8. continued. 



Element 
(Ionization) 


AiA) 




Accuracy 


references 


Element 
(Ionization) 


a 

AiA) 




Accuracy 


references 


Nil 


4604.9820 


-0.250 


D 


Kurucz 


Srll 


4077.7090 


0.151 




NIST 


Nil 


4606.2190 


-1.000 


D 


FMW 


Srll 


4215.5200 


-0.169 




NIST 


Nil 


4609.9050 


-0.580 




Kurucz 












Nil 


4617.8620 


-0.525 




Kurucz 


YII 


4883.6840 


0.070 




kurucz 


Nil 


4631.0170 


-0.957 




Kurucz 


YII 


4900.1200 


-0.09 




kurucz 


Nil 


4648.6460 


-0.100 


D 


Kurucz 


YII 


4982.1290 


-1.290 




kurucz 


Nil 


4886.7050 


-1.780 




Kurucz 


YII 


5087.4160 


-0.170 




kurucz 


Nil 


4886.9760 


-1.120 




Kurucz 


YII 


5200.4060 


-0.570 




kurucz 


Nil 


4900.9670 


-1.670 


E 


FMW 












Nil 


4904.4070 


-0.170 


D 


FMW 


Zrll 


4149.2170 


-0.030 




kurucz 


Nil 


4912.0200 


-0.800 


D 


FMW 


Zrll 


4156.2400 


-0.776 




kurucz 


Nil 


4913.9680 


-0.630 


D 


FMW 


Zrll 


4161.2100 


-0.720 




kurucz 


Nil 


4918.3620 


-0.240 


D 


FMW 


Zrll 


4208.9800 


-0.460 




kurucz 


Nil 


4918.7060 


-0.780 




Kurucz 


Zrll 


4496.9600 


-0.810 




Bal81 


Nil 


4925.5590 


-0.770 


D 


Kurucz 












Nil 


4935.8310 


-0.350 


D 


FMW 


Ball 


4554,0290 


0,163 


B 


Kurucz 


Nil 


4937.3410 


-0.390 


D 


FMW 


Ball 


4934,0760 


-0,156 


B 


NBS 


NU 


4976.3260 


-3.100 


C-H 


NIST 


Ball 


5853,6680 


-1,510 


B 


NIST 


Nil 


5003.7410 


-2.800 


C-i- 


NIST 


Ball 


6141,7130 


-0,810 


B 


NIST 


Nil 


5080.5280 


0.330 




kurucz 


Ball 


6496,8970 


-1,010 


B 


NIST 


Nil 


5081.1070 


0.300 




kurucz 












Nil 


5084.0890 


0.090 




kurucz 












Nil 


5096.8540 


-0.900 




kurucz 












Nil 


5099.9270 


-0.100 




kurucz 












Nil 


5102.9660 


-2.620 


c+ 


NIST 












Nil 


5137.0740 


-1.990 


c+ 


NIST 












Nil 


5476.9040 


-0.890 


c+ 


NIST 












Nil 


5553.6900 


-3.240 


c+ 


NIST 












Nil 


5587.8580 


-2.140 


C-i- 


NIST 












Nil 


5592.2620 


-2.590 


c+ 


NIST 












Nil 


5711.8880 


-2.260 


c+ 


NIST 












IN 11 


J / J4.0J0U 


'3/1 r\ 




iNlo 1 












Nil 


5892.8720 


-2.340 


C-H 


NIST 












Nil 


6163.4180 


-0.682 




kurucz 












Nil 


6170.5670 


-1.808 




kurucz 












Nil 


6175.3600 


-0.530 




kurucz 












Nil 


6176.8070 


-0.260 




kurucz 













